Supernova light-curve fitters and Dark Energy 
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We show that when a procedure is made to remove the tension between a supernova la (SN la) 
data set and observations from BAO and CMB, there might be the case where the same SN la 
set built with two different light-curve fitters behaves as two separate and distinct supernova sets, 
and the tension found by some authors between supernova sets actually could be due to tension or 
inconsistency between fitters. We also show that the information of the fitter used in an SN la data 
set could be relevant to determine whether phantom type models are favored or not when such a 
set is combined with the BAO/CMB joint parameter. 



I. INTRODUCTION 

The type la supernova (SN la) measurements remain 
a key ingredient in all current determinations of cosmo- 
logical parameters. More than a decade ago, combined 
observations of nearby and distant SNe la led to the dis- 
covery of the accelerating universe picture. It has be- 
come clear that different cosmological observations, such 
as the dimming of distant SNe la pHlj], anisotropics in 
the cosmic microwave background (CMB) [J .and the 
signature of baryon acoustic oscillation (BAO) HIB] can- 
not be explained with a cosmological model that contains 
only baryonic and dark matter according to a FRW stan- 
dard model. The most popular solution is to introduce 
an extra component with negative pressure, the so-called 
dark energy (e.g., @-[llj|). 

Characterization of dark energy focuses on estimation 
of the equation of state w, which is the ratio of pressure 
to density. For the time- invariant w — — 1, the equation 
of state is consistent with a cosmological constant. Any 
other fixed or time variable value of w would require more 
exotic models. 

It is a known fact that the same SN la data set in which 
distance estimates are analyzed with two different light- 
curve fitters, the values achieved for various cosmological 
parameters (for example, the equation of state of dark 
energy) differ, or also there could be found that some 
cosmological models result more favored than others (e.g. 

EM3)- 

Here we analyze the consistency between the fitters 
MLCS2k2 pj] (hereafter MLCS) and SALT2 p]| from a 
particular approach as it will be further explained. The 
Multicolor Light Curve Shape fitter, MLCS, is the most 
recent incarnation of the fitter used by the High-z Super- 
nova Team [2] , whilst the Spectral Adaptive Light curve 
Template, SALT2, is an improved version of the fitter 
used originally by the Supernova Cosmology Project p|. 
A detailed description of both fitters and a thorough dis- 
cussion about systematic errors in SN surveys can be 
found for example in pj|, . 
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Each method results in a distance modulus for each 
supernova. However, distance moduli calculated for the 
same objects by the two fitting methods are not necessar- 
ily equal. Whereas the MLCS calibration uses a nearby 
training set of SNe la assuming a close to linear Hubble 
law, SALT2 uses the whole data set to calibrate empirical 
light curve parameters. SNe la from beyond the range in 
which the Hubble law is linear are used, so a cosmologi- 
cal model must be assumed in this method. Typically a 
ACDM or a wCDM (w — const) model is assumed. Con- 
sequently the published values of SN la distance moduli 
obtained with SALT2 fitter retain a degree of model de- 
pendence. Regarding this, in |20j it was pointed out that 
systematic errors in the method of SNe la distance esti- 
mation have come into sharper focus as a limiting factor 
in SN cosmology. The major systematic concerns for su- 
pernova distance measurements are errors in correcting 
for host-galaxy extinction and uncertainties in the intrin- 
sic colors of supernovae, luminosity evolution, and selec- 
tion bias in the low-redshift sample. Also, SALT2 fitter 
does not provide a cosmology-independent distance esti- 
mate for each supernova, since some parameters in the 
calibration process are determined in a simultaneous fit 
with cosmological parameters to the Hubble diagram. 

In [2l[ it was investigated the tension between SN la 
data sets (Constitution pjj and Union (22[) and other 
data sets, CMB and BAO. There, it was shown that SN la 
data sets are in tension not only with the observations of 
CMB and BAO, but also with other SN la data sets such 
as Davis07 [23j . It was also shown that in the Davis07 
data set there is no tension with CMB and BAO observa- 
tions, concluding that Union and Constitution data sets 
are in tension not only with the observations of CMB and 
BAO, but also with other SN la data sets. Then, the au- 
thor found the main sources responsible for the tension 
by employing a truncation method, following the simple 
procedure used in [24f . With this in mind, the truncated 
UnionT and ConstitutionT data sets were built, which 
are consistent with the other observations and the ten- 
sion was completely removed. 

In the quest of the characterization of dark energy, in 
[25j it was first noted that observational data do not rule 
out the possibility that w < — 1. Phantom dark energy 
models with w < — 1 have the interesting properties that 



2 



the density of the dark energy increases with increasing 
scale factor, and the phantom energy density can become 
infinite at a finite time, a condition known as the 'big rip'. 

Future projects as SNAP (2(| are designed to reveal the 
nature of the dark energy. It will characterize the dark 
energy density, equation of state and time variation by 
precisely and accurately measuring the distance-redshift 
relation of SNe la. The matter density, dark energy den- 
sity, and flatness of the universe could be determined at 
the 1% level, including systematic uncertainties, the dark 
energy equation of state to about 3% and its time varia- 
tion characterized to within 10% of the Hubble expansion 
time. In summary, this kind of project will seek to de- 
termine whether w = — 1 or not, and if w is a constant 
or it evolves in time. However, in order to do this, the 
matter of reducing the systematic uncertainties between 
light-curve fitters will be of great importance. The com- 
bination of observational data sets might lead, as it will 
be further shown, the equation of state of dark energy 
to be of the phantom type or not, depending on how the 
SN la data were processed. 

In this Letter we focus on the consistency of two of 
the main light-curve fitters used for the elaboration of 
SN la data sets. To accomplish this, we present an- 
other approach to the tension between SN la data sets 
as the one found in 2l| and then the consistency be- 
tween light-curve fitters is analyzed when appears the 
need of removing tension between an SN la data set and 
BAO/CMB data applying a truncation method to the 
same SN la data set, but obtained with the two different 
fitters. Additionally, we show that the conclusion about 
if the combination of SN la data with BAO/CMB favors 
or not an equation of state for dark energy of the phan- 
tom type in the framework of a given cosmological model 
might depend, in some cases, on the fitter employed in 
the elaboration of the SN la set used in the analysis. 



II. TENSION BETWEEN SUPERNOVAE IA 
DATA SETS REVISITED 

In this section and the next one we will use a x 2 = 
XsNe+XsAO/CMB statistic to analyze the confidence in- 
tervals of the free parameters of two cosmological models, 
by employing different SN la data sets and their combi- 
nation with the BAO/CMB joint parameter introduced 
in 14]. Since lately more non-standard models are built 
to try to explain the dark energy phenomenon, we chose 
to combine the SNe analysis with the BAO/CMB joint 
parameter which is better suited for these sort of mod- 
els, using it only as an example of data combination. We 
must also clarify that the goal of this work was neither to 
put constraints nor to find best fits to cosmological mod- 
els, but to show certain extra information that should be 
minded when non-SNe data sets are combined with SNe 
data sets analyzed with different fitters. The separate \ 2 
of SNe la and BAO/CMB used in this work are shown 
in Appendix A. 



We considered for the analysis the ACDM and flat 
wCDM models with w — const. In the case of ACDM the 
dark energy is a cosmological constant which behaves as 
a vacuum energy with w = — 1, but we allowed non-zero 
spatial curvature flk- Then the dimensionless expansion 
rate is given by, 

E ee H(z)/H Q = 

[o m (i + z) 3 + o r (i + zf + n fc (i + zf + n A } 1/2 (i) 

where H(z) is the Hubble parameter as a function of the 
cosmological redshift z; VL m , Q r and £Ia are the contribu- 
tions of matter, radiation and dark energy respectively 
to the total energy density today, and the curvature den- 
sity is fife = 1 — £l m — f2 r — J7a- The parameters usually 
chosen as free parameters in this model are VL m and Q\. 

For the flat wCDM (fi^ = 0) case, we allowed the 
equation of state parameter of dark energy w to differ 
from -1 so, 

E = [n m (i + z) 3 + n r (i + z) 4 + n A {i + zf ( - 1+w '>} 1 / 2 (2) 

where now J7a = l — ^l m — ^r- The free parameters in this 
model are Q m and w. The case flat ACDM is a special 
case of this one. 

The choice of the considered models for this analysis 
was arbitrary and not relevant, this was for mere sim- 
plicity (two free parameters) with the aim of showing the 
novel results. The consideration of models with an equa- 
tion of state w that does not evolve in time is enough, 
since the current data do not yield precise constraints on 
the time derivative of w. (For constraints on time- varying 
w, see [HI as an example). Additionally, the rationale for 
the flat case of wCDM was that the WMAP data from 
the CMB anisotropy constrain the spatial curvature to be 
very small (e.g. [27|). Also, the use of these two models is 
the most suitable to obtain a correct comparison among 
different fitters, since the light-curve fitter SALT2 uses 
the whole data set to calibrate empirical light curve pa- 
rameters and a ACDM or a wCDM model is typically 
assumed, as it was mentioned in Section I. 

As it was previously mentioned, the analysis in this 
work was performed in the framework of SALT2 [l9| and 
MLCS [18[ fitters and the SN la data sets used were 
Constitution data sets (Tables 2 and 4 from [H] as it 
will be further explained in Section III), the SDSSII full 
data set (Tables 10 and 14 from [l3[ with the same values 
used for the 'intrinsic' dispersions there) and the Union2 
data set [|[ . 

It was interesting to observe the discrepancy between 
the results obtained when using one fitter or the other 
when it is allowed the variation of flk in the framework 
of a ACDM (w = -1) model. The analysis of the SDSSII 
full data set with 288 SNe la built with MLCS showed 
that the flat case (Cl m — 0.27 [H) stands excluded to 
more than 3 a confidence level, while with the same data 
set, but processed with SALT2 this did not happen at 
all (Figs, la and lc). Something similar happens in the 
framework of the flat wCDM model (Figs, lb and Id): 
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the standard model (fi m = 0.27, w — —1) using SDSSII 
(MLCS) is excluded to more than 2 a confidence level 
(see also Fig.l of [Hj]). Since the responsible of this fact 
is the fitter and not the SNe la (because the data set 
is the very same and only the fitter was changed), one 
could then wonder what SN la data set should be used to 
be combined with, for example CMB data, which leave 
little margin to the variation of (e.g. j27|). Looking 
at Figs, lc and Id one would choose those data sets 
processed with SALT2; however we should keep in mind 
that SALT2 fitter retains a degree of model dependence 
because typically a ACDM model is assumed. 
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FIG. 1: (a) Confidence intervals at 68.3%, 95.4% and 99.7% 
in the Q, m - fl k plane for the SDSSII (MLCS) SN la data 
set in the ACDM model framework (dashed lines), (b) Con- 
fidence intervals at 68.3% and 95.4% in the Q, m — w plane for 
the SDSSII (MLCS) SN la data set in the flat raCDM model 
framework (dashed lines), (c) Confidence intervals at 68.3% 
and 95.4% in the Q, m - Q k plane for the SDSSII (SALT2) SN 
la data set in the ACDM model framework (dashed lines), (d) 
Confidence intervals at 68.3% and 95.4% in the Q m —w plane 
for the SDSSII (SALT2) SN la data set in the flat wCDM 
model framework (dashed lines). The best fits are indicated 
with a star whereas the standard flat ACDM (Q m = 0.27) is 
marked with a triangle. 

In l2ll . it was found that Union and Constitu- 
tion [I2| data sets were in tension not only with BAO 
and CMB data, but with other SN la data sets too 
(e.g. Davis07 [23]). The tension found was attributed 
to certain supernovae of the data set and by a truncation 
method these outliers were removed from the set with 
the objective of releasing the tension. 

We found interesting to analyze what would happen if 
we applied a criterion in order to study the consistency 
between data sets, a criterion more restrictive than the 
only fact that the confidence intervals overlap. To per- 
form this analysis, we adopted the criterion of consider- 
ing the existence of tension between a given data set and 



another set constituted combining several data sets (in- 
cluding the first one) as the fact that the best fit point to 
the first data set is out of the 68.3% (la) confidence level 
contour given by the combined data set. Similar criteria 
were adopted in their analysis by [H [H [H HI] . This 
is a criterion we will adopt in order to show how differ- 
ently will behave the same SN la set processed with two 
different fitters when a truncation method is performed. 
This will lead to an alternative way of analyzing the dis- 
crepancy between the results obtained when one or other 
fitter is used. One could choose not to use this more 
restrictive criterion, nevertheless with this adopted crite- 
rion, we seek more physical consistency between best-fits, 
so the best fits do not drive to too different cosmological 
evolutions. A best fit which effective equation of state is 
of the phantom type [2|| (w < —I) tells us about very 
different physics from the one that is not. For instance, 
in a recent work [2^ | the consequences of applying it to 
several data sets in the framework of f(T) theories have 
been investigated. 

When we used the Union2 recently released data set 
with 557 SNe la (processed with SALT2 fitter) to com- 
bine it with the information from BAO/CMB in the 
framework of the flat wCDM model, we found that there 
is no tension between data sets (Fig. 2a). However, when 
we did the same procedure with the 288 SNe la of the 
SDSSII set (in the framework of MLCS fitter) wc found 
that there is tension between SNe la and BAO/CMB to 
more than 3 er (Fig. 2b). Since the data from BAO/CMB 
are the same used in both cases, one could wrongly con- 
clude that there is tension between the SN la data sets. 
This case is similar to what was found in 12 II . The result 
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FIG. 2: (a) Confidence intervals at 68.3% and 95.4% in 
the Q m - w plane for the Union2 (SALT2) SNe la only 
(dashed lines) and Union2+BAO/CMB (solid lines), (b) Con- 
fidence intervals at 68.3%, 95.4% and 99.7% in the Q m — w 
plane for the SDSSII (MLCS) SNe la only (dashed lines) and 
SDSSII+BAO/CMB (solid lines). The best fits to SNe are 
indicated with a star whereas the best fits to the combined 
parameters are indicated with a dot. 
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obtained in [21], where there is tension between Union 
and BAO, while in Davis07 there is none, seems to be 
due to a data set (Union) is processed with SALT and 
the other (Davis07) with versions of MLCS. The tension 
between Union and Davis07 could actually be a trouble 
between fitters and not between supernovae sets as it will 
be discussed in the next section. 

To better understand the analysis of tension between 
Union2 and SDSSII performed here, in Figs. 3a and 3b 
are displayed the combinations of SNe Ia+BAO/CMB 
for the Union2 and SDSSII cases respectively. For the 
case Union2 vs SDSSII(MLCS) there is tension to more 
than 2 a level, while in Union2 vs SDSSII(SALT2) there 
is none. Clearly the tension exists between fitters and 
not between SN la data sets. 
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FIG. 3: Confidence intervals at 68.3% and 95.4% in the 
fl m —w plane from combining SNe Ia+BAO/CMB. (a) Union2 
(SALT2 - dashed lines) vs SDSSII (MLCS - solid lines), (b) 
Union2 (SALT2- dashed lines) vs SDSSII (SALT2 - solid 
lines). A tension between fitters is revealed. 



III. ARE SUPERNOVAE LIGHT-CURVE 
FITTERS CONSISTENT? 

In this section we analyze the problem of tension be- 
tween light-curve fitters and the consistency between 
them from a different approach. 

We found interesting to study what happens when one 
needs to perform a truncation procedure as in [2l[ to 
remove the tension between BAO/CMB and an SN la 
data set (processed with a given fitter) and the same SN 
la data set but with a different fitter. It is important to 
emphasize that the adopted criterion itself is not relevant, 
but how the same SN set behaves when it is processed 
with different fitters. 

The criterion consists in finding and removing the out- 
liers responsible of the tension. First, we fitted the model 
to the whole SN la data set finding the best fit parame- 
ters, including the nuisance parameter (Iq (see Appendix 



A). Then, we calculated the relative deviation to the best 
fit prediction, | fj, f, s ~ n t h | /fobs; for all the data points 
finding which cut solved the tension problem and which 
SNe la were the outliers. 

We performed this procedure in the framework of the 
flat wCDM model for two different SN la data sets, Con- 
stitution (SALT2 and MLCS) and SDSSII (SALT2 and 
MLCS). 

For the data sets here denominated Constitution 
(SALT2) and Constitution (MLCS) we used the same 337 
SNe la from the Table 2 (SALT2) and Table 4 (MLCS17) 
from H3. In 112] MLCS was used to find the dust- 
reddening properties through the value of Rv that min- 
imizes the scatter in the Hubble residuals for the nearby 
CfA3 sample and they found R v = 1.7. MLCS with 
Rv — 3.1 overestimates host-galaxy extinction while 
R v = 1.7 does not. With its lower value of R v , MLCS17 
attributes less host extinction to each SN la and therefore 
this produces a larger distance compared to MLCS31. 

Figure 4a shows the confidence intervals to 68.3% and 
95.4% in the plane fi TO - w for the Constitution (SALT2) 
data set of SNe la only and for the combination SNe 
Ia+BAO/CMB. There, a tension to 2 a confidence level 
between both data sets can be observed. In Fig. 4b it 
is shown that after the truncation procedure with a 2.1 
a cut (13 outliers) the tension was removed completely. 
In a similar way, we proceeded to do the same with the 
Constitution (MLCS) data set and the results are shown 
in Figs. 5a and 5b. In this case, a 2 a cut was enough and 
12 SNe la were removed from the data set so the tension 
with BAO/CMB was removed. Although the necessary 
cut for Constitution set was slightly different for SALT2 
than for MLCS, and although the number of SNe to re- 
move was not significantly different, the most remarkable 
fact was that only 4 outliers were the same. This can be 
clearly seen in Table[T]where SNe la outliers are displayed 
for each case and which of those are the ones in common. 
We want to stress that the SN la data set was the same 
and the only difference lied in the light-curve fitter em- 
ployed. Then, with this truncation procedure we found 
the behavior was as if there were two different SN la sets 
when actually there was only one. 

Something more drastic occurred when we did the 
same analysis, but using the SDSSII data set. When the 
SDSSII (SALT2) data set was used, no SNe la needed to 
be removed, since there was no tension with BAO/CMB 
(Fig. 6a). However, with the SDSSII (MLCS) data set 
there was tension with BAO / CMB to a level greater than 
99.7% (Fig. 6b) and a cut of 1.7 a (18 outliers) was 
needed to remove such tension. To highlight this, in Fig. 
7 we show the confidence intervals to 99.7% of the com- 
bined SNe Ia+BAO/CMB for the cases SDSSII (SALT2) 
and SDSSII (MLCS). There, it can be seen how both 
best fits differ by more than 3 a level. Another inter- 
esting thing is the comparison of Figs 7 and 3a. This 
comparison allows appreciating how two light-curve fit- 
ters employed for the same SN la set produce the same 
result than two different SN la sets. 



5 



-0.5- 
-LO- 



CI. 3 

n 



FIG. 4: (a) Confidence intervals at 68.3% and 95.4% in the 
f2 m — w plane for the Constitution (SALT2) SNe la only 
(dashed lines) and SNe+BAO/CMB (solid lines). Tension 
can be appreciated between data sets, (b) Idem (a), but after 
performing the truncation procedure. The best fits to SNe are 
indicated with a star whereas the best fits to the combined 
parameters are indicated with a dot. 
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5: (a) Confidence intervals at 68.3% and 95.4% in the 
w plane for the Constitution (MLCS) SNe la only 
(dashed lines) and SNe+BAO/CMB (solid lines). Tension 
can be appreciated between data sets, (b) Idem (a), but after 
performing the truncation procedure. The best fits to SNe are 
indicated with a star whereas the best fits to the combined 
parameters are indicated with a dot. 



The fact of needing to take out different supernovae 
to remove the tension in Constitution (SALT2) and in 
Constitution (MLCS) reveals that they are not the su- 
pernovae themselves that probably generate the tension, 
but the fitters. Using the same SN la data set with the 
same truncation process, the SNe la we had to take out 
were not the same ones. It seems as if a given fitter makes 
some supernovae bring an apparent problem. The case 



TABLE I: The outliers after the truncation procedure by 
using the same Constitution SN la data set processed with 
MLCS and SALT2 light-curve fitters. Note that there are 
only 4 matches in the outlier SNe. 



SALT2 



04D3cp 



£221 


£221 


isis 


isis 



sn06cg 



sn07bz 



silMhr 



sn97dg 
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FIG. 6: (a) Confidence intervals at 68.3% and 95.4% in the 
0„ - w plane for the SDSSII (SALT2) SNe la only (dashed 
lines) and SNe+BAO/CMB (solid lines). There is no tension 
between the data sets, (b) Idem (a), but for the SDSSII 
(MLCS) SNe la. Here can be appreciated tension to more 
than 2 a confidence level. 



of the SDSS data set also presented this behavior, be- 
cause with one fitter there were supernovae that caused 
tension and with the other there were none, since there 
was no need of truncation. What we found here with two 
different fitters for the same SN la data set is analogue 
to what was found in [2l| between two different SN la 
sets. What was found between Union and Davis07 data 
sets is perhaps a fitter problem and not supernovae. 

We also analyzed another interesting aspect. We found 
important to stress how an employed fitter in an SN la 
data set could take part in the conclusion about the 
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FIG. 7: Confidence intervals at 68.3% and 99.7% in the 
fl m — w plane coming from combining SN la and B AO / CMB 
data. Solid lines correspond to SDSSII (SALT2) whereas 
dash-dotted lines correspond to SDSSII (MLCS). There can 
be seen that both best fits are outside 3 a confidence level 
from the other fitter. 



equation of state w of a given model when this SN 
la set is combined with BAO/CMB, being this joint 
parameter more suitable for the lately more common 
non-standard models. In Table |TT] (Phantom Table) we 
show how for two given theories using the combination 
BAO/CMB+SN la (SALT2) is obtained a best fit for w 
of the phantom type, while with MLCS the opposite oc- 
curs. There, the result is shown for the flat wCDM model 
and for a non-standard model of modified gravity f(T) as 
in [H,[30| when Union2 (SALT2), SDSSII (SALT2 and 
MLCS) and Constitution (SALT2 and MLCS) are used. 



TABLE II: Phantom Table. Best fit w types for two cosmolog- 
ical models when the BAO/CMB joint parameter is combined 
with an SN la data set processed with a given light-curve fit- 
ter. 



Model 


BAO/CMB + 


Best fit it; type 


Flat wCDM 


Union2 (SALT2) 


Phantom 


Flat wCDM 


SDSSII (SALT2) 


Phantom 


Flat wCDM 


SDSSII (MLCS) 


non-Phantom 


Flat wCDM 


Const.(SALT2) 


Phantom 


Flat wCDM 


Const. (MLCS) 


non-Phantom 


f(T) 


Union2 (SALT2) 


Phantom 


}{T) 


SDSSII (SALT2) 


Phantom 


f(T) 


SDSSII (MLCS) 


non-Phantom 


f(T) 


Const.(SALT2) 


Phantom 


f(T) 


Const. (MLCS) 


non- Phantom 



We do not mean that this is going to occur with any 
model, but to conclude that the combination of data sets 
favors or not phantom type models with w < — 1, the fit- 
ter used to process the SNe la is an additional factor that 
must be taken into account as source of degeneration. 



IV. CONCLUSIONS 

Although the evidence in the FRW framework that 
the universe is going through an accelerated stage, be- 
cause of the existence of what we call dark energy, is 
solid from various observational data sets, its nature will 
depend on the future better understanding of the sys- 
tematic errors present in supernovae observations, par- 
ticularly those present in the light-curve fitters analyzed 
here. 

As it is well known, when the same SN la data set is 
processed with two different fitters, the values found for 
cosmological parameters (such as the equation of state of 
dark energy) differ. 

Here we analyzed this difference showing how when a 
procedure to remove tension between an SN la data set 
and observations from BAO/CMB is performed, there 
could exist the case where the same SN la set processed 
with two different fitters behaves as if there were two dif- 
ferent sets, and the tension between sets found by some 
authors actually could be due to a tension or inconsis- 
tency between fitters. 

We also showed that the information of the fitter used 
in an SN la set could be relevant and it should be minded 
as an additional factor to decide if phantom type models 
are favored or not when the given SN la set is combined 
with the BAO/CMB joint parameter. 
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Appendix A: Cosmological constraints methods 

1. Type la Supernovae constraints 

The N data points of the SNe la compiled in a data 
set are usually given in terms of the distance modulus 
l^obs(zi)- On the other hand, the theoretical distance 
modulus is defined as 



Hth{zi) = blog w D L (zi) + /i 



(Al) 



where no = 42.38 — 5log\oh and h is the Hubble constant 
Hq in units of 100 km/s/Mpc, whereas the Hubble-free 
luminosity distance for the general case is, 



# L (z) = (i + z)|ft fc r 1/2 sJ|r! fc | 1/2 



dz' 



(A2) 



in which E = H(z)/Hq is the dimensionless expansion 
rate, p denotes the model parameters, and the function 



7 



Sk(x) = sin(x) when the curvature density ttk < 0, 
<Sfc(#) = sinh(x) for J\. > and Sk(x) = x for the flat 
case fife = 0. Correspondingly, the \ 2 from the N SNe la 
is given by 



XSNeKP) — /^i -K „ \ \ A6 > 



where a(zi) is the corresponding uncertainty for each ob- 
served value. The parameter fj,Q is a nuisance parameter 
but it is independent of the data points. One can per- 
form an uniform marginalization over (1q. However, there 
is an alternative way. Following [31], the minimization 
with respect to /xq can be made by expanding the XsNe 
of (| A3|) with respect to fj,Q as 



XsNe(p) =A-2imiB + i4C 



(A4) 



where, 



N 



Mp) = E 



i=l 
N 



B(P) = E 



[Hobs(zj) - Hth(Zi-,Ho = 0, p)] 2 

[Hobs(zi) ~ nth(zi;no = 0,p)] 



i=l 
2V 



CJ 2 {Zi) 



C = T — 



Eq. (|A4j) has a minimum for /Xo = S jC at 

B(p) 2 



XSiVe(P) = MP) ~ 



c 



(A5) 



Since XsNe, min = XsNe,min obviously, we can instead 
minimize XsNe which is independent of hq. 

2. Combined BAO/CMB parameter constraints 

Since in this work we combined CMB and BAO obser- 
vations with SN la data sets for flat models, we consider 
here all the relations for the spatially flat case. 

A more model-independent constraint can be achieved 
by multiplying the BAO measurement of r s (z^)j Dy(z) 
with the position of the first CMB power spectrum peak 
27] I a = ^d,A{z Sf )/r s (z^), thus canceling some of the de- 
pendence on the sound horizon scale [l4|. Here, cLa(z*) 
is the comoving angular-diameter distance to recombi- 
nation, r s is the comoving sound horizon at photon de- 
coupling, Zd w 1020 is the redshift of the drag epoch at 



which the acoustic oscillations are frozen in, and Dy is 
defined as (assumed a ACDM model) @, 



D v (z) 



C/n H(z')) 



1 1/3 



(A6) 



H(z) \J Q H(z') 
We further assume z* = 1090 from [27] (variations within 
the uncertainties about this value do not give significant 
differences in the results). 

In @ was measured r s {zd) / 'Dy(z) at two redshifts, 
z = 0.2 and z = 0.35, finding r s (z d )/ D v (0.2) = 0.1905 ± 
0.0061 and r s (z d )/D v (0.35) = 0.1097 ± 0.0036. Com- 
bining this with I a gives the combined BAO/CMB con- 
straints flij : 



d A {z*) r s {z d ) 



D v (0.2) r s (z*) 
d A {z*) r s (z d ) 



LV(0-35) r s (z* 



18.32 ±0.59 
10.55 ±0.35 



(A7) 



Before matching to cosmological models we also need 
to implement the correction for the difference between 
the sound horizon at the end of the drag epoch and the 
sound horizon at last scattering. The first is relevant for 
the BAO, the second for the CMB, and r s (z d )/r s (z*) — 
1.044 ±0.019 (using values from (27]]). Inserting this into 
(|A7|) and taking into account the correlation between 
these measurements using the correlation coefficient of 
0.337 calculated by @, gives the final constraints we use 
for the cosmology analysis [lij ]: 



At 
A 2 = 



<Ia{z*) 
= Dv(0.2) 

<Ia{z*) 
LV(0.35) 



17.55 ±0.65 

10.10 ±0.38 (A8) 



Using this BAO/CMB parameter cancels out some of 
the dependence on the sound horizon size at last scat- 
tering. This thereby removes the dependence on much 
of the complex pre-recombination physics that is needed 
to determine that horizon scale [l4|. In all the cases, we 
have considered a radiation component f2 r =5 x 10 . 

So, for our analysis we add to the % 2 statistic: 



Xbao/Cmb (P) 



N = 2 r . , , 

y-> [Aob s (Zi) 



A th (zi;p)f 



o\(Zi) 



(A9) 



where p are the free parameters, A a b s is the observed 
value (Ai and A 2 ), A th is the predicted value by the 
model and a a is the la error of each measurement. 
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